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Abstract 

d ' The so-called Klein paradox - unimpeded penetration of relativistic particles 

J^ . through high and ^vide potential barriers - is one of the most exotic and counterintu- 

^. . . . 

itive consequences of quantum electrodynamics (QED). The phenomenon is discussed 

"t^ I in many contexts in particle, nuclear and astro- physics but direct tests of the Klein 

paradox using elementary particles have so far proved impossible. Here w^e show 

^ . that the effect can be tested in a conceptually simple condensed-matter experiment 

O . 

CJ \ by using electrostatic barriers in single- and bi-layer graphene. Due to the chiral 

—^ ■ nature of their quasiparticles, quantum tunneling in these materials becomes highly 

>■..... 
^r^ • anisotropic, qualitatively different from the case of normal, nonrelativistic electrons. 

<N , 

ff^ , Massless Dirac fermions in graphene allow a close realization of Klein's gedanken 

^^ ' experiment w^hereas massive chiral fermions in bilayer graphene ofTer an interesting 

^: 

^^ ■ complementary system that elucidates the basic physics involved. 
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The term Klein parado x^i^i^i"^i^i^i^ usually refers to a counterintuitive relativistic process 
in which an incoming electron starts penetrating through a potential barrier if its height 
Vq exceeds twice the electron's rest energy mc^ (where m is the electron mass and c the 
speed of light). In this case, the transmission probability T depends only weakly on the 
barrier height, approaching the perfect transparency for very high barriers, in stark contrast 
to the conventional, nonrelativistic tunneling where T exponentially decays with increasing 
Vq. This relativistic effect can be attributed to the fact that a sufficiently strong potential, 
being repulsive for electrons, is attractive for positrons and results in positron states inside 
the barrier, which align in energy with the electron continuum outside^>^. Matching be- 
tween electron and positron wavefunctions across the barrier leads to the high-probability 
tunneling described by the Klein paradox^. The essential feature of QED responsible for 
the effect is the fact that states at positive and negative energies (electrons and positrons) 
are intimately linked (conjugated), being described by different components of the same 
spinor wavefunction. This fundamental property of the Dirac equation is often referred to 
as the charge- conjugation symmetry. Although Klein's gedanken experiment is now well 
understood, the notion of paradox is still used widely^'^'^'^^, perhaps because the effect has 
never been observed experimentally. Indeed, its observation requires a potential drop ^ mc^ 
over the Compton length h/mc, which yields enormous electric f\.e\d^^{S > lO^^V/cm) and 
makes the effect relevant only for such exotic situations as, for example, positron production 
around super- heavy nuclei^ with charge Z > 170 or evaporation of black holes through 
generation of particle-antiparticle pairs near the event horizon^. The purpose of this pa- 
per is to show that graphene - a recently found allotrope of carbon^ - provides an effective 
medium ("vacuum") where relativistic quantum tunneling described by the Klein paradox 
and other relevant QED phenomena could be tested experimentally. 

DIRAC-LIKE QUASIPARTICLES IN GRAPHENE 

Graphene is a single layer of carbon atoms densely packed in a honeycomb lattice, or it 
can be viewed as an individual atomic plane pulled out of bulk graphite. From the point 
of view of its electronic properties, graphene is a two-dimensional zero-gap semiconductor 
with the energy spectrum shown in Fig. la and its low-energy quasiparticles are formally 
described by the Dirac-like HamiltonianiSiiLi^ 

Hq = —ihvFcrV (1) 



where vp ~ 10^ ms~^ is the Fermi velocity, and a = (cr^., ay) are the Pauh matrices. Ne- 
glecting many-body effects, this description is accurate theoreticalljsiSiiiJ^ and has also been 
proven experiment allyi^i^^ by measuring the energy- dependent cyclotron mass in graphene 
(which yields its linear energy spectrum) and, most clearly, by the observation of a relativistic 
analogue of the integer quantum Hall effect. 

The fact that charge carriers in graphene are described by the Dirac-like equation (1) 
rather than the usual Schrodinger equation can be seen as a consequence of graphene's crystal 
structure, which consists of two equivalent carbon sublattices A and B^°i^^i^^ . Quantum 
mechanical hopping between the sublattices leads to the formation of two cosine-like energy 
bands, and their intersection near the edges of the Brillouin zone (shown in red and green in 
Fig. ^) yields the conical energy spectrum. As a result, quasiparticles in graphene exhibit 
the linear dispersion relation E = hkvp, as if they were massless relativistic particles (for 
example, photons) but the role of the speed of light is played here by the Fermi velocity 
vp ^ c/300. Due to the hnear spectrum, one can expect that graphene's quasiparticles 
behave differently from those in conventional metals and semiconductors where the energy 
spectrum can be approximated by a parabolic (free-electron-like) dispersion relation. 

Although the linear spectrum is important, it is not the only essential feature that under- 
pins the description of quantum transport in graphene by the Dirac equation. Above zero 
energy, the current carrying states in graphene are, as usual, electron-like and negatively 
charged. At negative energies, if the valence band is not full, its unoccupied electronic states 
behave as positively charged quasiparticles (holes), which are often viewed as a condensed- 
matter equivalent of positrons. Note however that electrons and holes in condensed matter 
physics are normally described by separate Schrodinger equations, which are not in any way 
connected (as a consequence of the Seitz sum rule^^, the equations should also involve differ- 
ent effective masses). In contrast, electron and hole states in graphene are interconnected, 
exhibiting properties analogous to the charge-conjugation symmetry in QED ^^i^^'^^ . For the 
case of graphene, the latter symmetry is a consequence of its crystal symmetry because 
graphene's quasiparticles have to be described by two-component wavef unctions, which is 
needed to define relative contributions of sublattices A and B in quasiparticles' make-up. 
The two-component description for graphene is very similar to the one by spinor wavefunc- 
tions in QED but the 'spin' index for graphene indicates sublattices rather than the real 
spin of electrons and is usually referred to as pseudospin a. 





FIG. 1: Tunneling through a potential barrier in graphene. (a)- Schematic diagrams of the spec- 
trum of quasiparticles in single-layer graphene. The spectrum is linear at low Fermi energies (<1 
eV). The red and green curves emphasize the origin of the linear spectrum, which is the crossing 
between the energy bands associated with crystal sublattices A and B. The three diagrams illus- 
trate schematically the positions of the Fermi energy E across the potential barrier of height Vq and 
width D shown in (b). The Fermi level (dotted lines) lies in the conduction band outside the barrier 
and the valence band inside it. The blue filling indicates occupied states. The pseudospin denoted 
by vector a is parallel (antiparallel) to the direction of motion of electrons (holes), which also 
means that a keeps a fixed direction along the red and green branches of the electronic spectrum. 
(c) - Low-energy spectrum for quasiparticles in bilayer graphene. The spectrum is isotropic and, 
despite its parabolicity, also originates from the intersection of energy bands formed by equivalent 
sublattices, which ensures charge conjugation, similar to the case of single-layer graphene. 



There are further analogies with QED. The conical spectrum of graphene is the result 
of intersection of the energy bands originating from sublattices A and B (see Fig. la) and, 
accordingly, an electron with energy E propagating in the positive direction originates from 



the same branch of the electronic spectrum (shown in red) as the hole with energy —E 
propagating in the opposite direction. This yields that electrons and holes belonging to 
the same branch have pseudospin a pointing in the same direction, which is parallel to 
the momentum for electrons and antiparallel for holes (see Fig. la). This allows one to 
introduce chiralityi^i, that is formally a projection of pseudospin on the direction of motion, 
which is positive and negative for electrons and holes, respectively. The term chirality is 
often used to refer to the additional built-in symmetry between electron and hole parts 
of graphene's spectrum (as indicated by color in Fig. 1) and is analogous (although not 
completely identical^i^) to the chirality in three-dimensional QED. 

KLEIN PARADOX REFORMULATED FOR SINGLE-LAYER GRAPHENE 
Because quasiparticles in graphene accurately mimic Dirac fermions in QED, this con- 
densed matter system makes it possible to set up a tunneling experiment similar to that 
analyzed by Klein. The general scheme of such an experiment is shown in Fig. 1, where we 
consider the potential barrier that has a rectangular shape and is infinite along the y-axis: 

fW = (''°'""^"''' (2) 

I otherwise. 

This local potential barrier inverts charge carriers underneath it, creating holes playing the 
role of positrons, or vice versa. For simplicity, we assume in (2) infinitely sharp edges, 
which allows a direct link to the case usually considered in QEEU'^i^'^iL^. The sharp-edge 
assumption is justified if the Fermi wavelength A of quasiparticles is much larger than the 
characteristic width of the edge smearing, which in turn should be larger than the lattice 
constant (to disallow Umklapp scattering between different valleys in graphene)^. Such a 
barrier can be created by the electric field effect using a thin insulator or by local chemical 
dopingSii^iii. Importantly, Dirac fermions in graphene are massless and, therefore, there 
is no formal theoretical requirement for the minimal electric field S to form positron-like 
states under the barrier. To create a well-defined barrier in realistic graphene samples with 
a disorder, fields S ~ 10^ V/ cm routinely used in experimental*^ should be sufficient, which 
is eleven orders of magnitude lower than the fields necessary for the observation of the Klein 
paradox for elementary particles. 

It is straightforward to solve the tunneling problem sketched in Fig. lb. We assume that 
the incident electron wave propagates at an angle (p with respect to the x axis and then try 



the components of the Dirac spinor ipi and ip2 for the Hamiltonian H = Hq + V (x) in the 
foUowing form: 



ipi{x,y) 



ip2 {x, y) 



^^ik^x _|_ ^g jfc^x^ ^tkyy ^ X < 0, 

Q ( piKxX-\-i(p rpp trCxX i<p\ ptkyy T <^ C] 

s' (ae*«-^+^^ - fe-^"-^-^^) e'''yy, 0<x < D, (3) 

g^^ikxX+ikyy+icj) X > D 



where kp = 27r/A is the Fermi wavevector, k^ = kp cos0 and ky = kp sincj) are the wavevec- 



tor components outside the barrier, q^ = \ {E — Vq) /h?Vp — k"^, 9 = tan"-*^ (ky/q^) is the 
refraction angle, s = signE, s' = sign {E — Vq). Requiring the continuity of the wavefunc- 
tion by matching up coefficients a, b, t, r, we find the foUowing expression for the reflection 
coefficient r 

o • jrf. • / T^N sin (f)- ss' sin 6 , ^ . 

r = 2?p *^sin n D) 1 4) 

^^" ^ ss' [e-^i^^ cos {(t) + 9) + ei^^ cos (0 - 9)] - 2i sin {q,D) ' ^ ' 

Fig. 12^1' shows examples of the angular dependence of transmission probability T = \t\ = 
1 — \r\ calculated using the above expression. In the limit of high barriers |Vo| ^ \E\, the 
expression for T can be simplified to 

T = ""^''^ (5) 

1 — cos^ (QxD) sin^ (J) 

Equations (4,5) yield that under resonance conditions q^D = ttN , A^ = 0, ±1,... the 
barrier becomes transparent (T = 1). More significantly, however, the barrier remains 
always perfectly transparent for angles close to the normal incidence = 0. The latter is 
the feature unique to massless Dirac fermions and directly related to the Klein paradox in 
QED. One can understand this perfect tunneling in terms of the conservation of pseudospin. 
Indeed, in the absence of pseudospin-fiip processes (such processes are rare as they require a 
short-range potential, which would act differently on A and B sites of the graphene lattice), 
an electron moving to the right can be scattered only to a right-moving electron state or 
left-moving hole state. This is illustrated in Fig. la, where charge carriers from the "red" 
branch of the band diagram can be scattered into states within the same "red" branch but 
cannot be transformed into any state on the "green" branch. The latter scattering event 




FIG. 2: Klein-like quantum tunneling in graphene systems. Transmission probability T through a 
100-nm-wide barrier as a function of the incident angle for (a) single- and (b) bi-layer graphene. 
The electron concentration n outside the barrier is chosen 0.5 x 10^^ cm~^ for all cases. Inside the 
barrier, hole concentrations p are 1 x 10^^ and 3 x 10^^ cm~^ for red and blue curves, respectively 
(such concentrations are most typical in experiments with graphene). This corresponds to the Fermi 
energy E of incident electrons ~ 80 and 17 meV for single- and bi-layer graphene, respectively, and 
A ~ 50 nm. The barrier heights Vq are (a) 200 and (b) 50 meV (red curves) and (a) 285 and (b) 
100 meV (blue curves). 



would require the pseudospin to be flipped. The matching between directions of pseudospin 
0" for quasiparticles inside and outside the barrier results in perfect tunneling. In the strictly 
one-dimensional case, such perfect transmission of Dirac fermions has been discussed in the 
context of electron transport in carbon nanotubes^^^ (see also ref. [19]). Our analysis 
extends this tunneling problem to the two-dimensional (2D) case of graphene. 

CHIRAL TUNNELING IN BILAYER GRAPHENE 

To elucidate which features of the anomalous tunneling in graphene are related to the 
linear dispersion and which to the pseudospin and chirality of the Dirac spectrum, it is in- 
structive to consider the same problem for bilayer graphene. There are both differences and 
similarities between the two graphene systems. Indeed, charge carriers in bilayer graphene 
have parabolic energy spectrum as shown in Fig. Ic, which means they are massive quasi- 
particles with a finite density of states at zero energy, similar to conventional nonrelativistic 
electrons. On the other hand, these quasiparticles are also chiral and described by spinor 
wavefunctions22.i^, similar to relativistic particles or quasiparticles in single-layer graphene. 



Again, the origin of the unusual energy spectrum can be traced to the crystal lattice of 
bilayer graphene with four equivalent sublattices^i. Although "massive chiral fermions" do 
not exist in the field theory their existence in the condensed matter physics (confirmed 
experimentally?^) offers a unique opportunity to clarify the importance of chirality in the 
relativistic tunneling problem described by the Klein paradox. In addition, the relevant 
QED-like effects appear to be more pronounced in bilayer graphene and easier to test ex- 
perimentally, as discussed below. 

Charge carriers in bilayer graphene are described by an off-diagonal Hamiltonian22i2i 

H„^-!^( ° ^'--"A (6) 

which yields a gapless semiconductor with chiral electrons and holes having a finite mass 
m. An important formal difference between the tunneling problems for single- and bi- 
layer graphene is that in the latter case there are four possible solutions for a given energy 
E = ±h'^kp/2m. Two of them correspond to propagating waves and the other two to 
evanescent ones. Accordingly, for constant potential Vi, eigenstates of Hamiltonian (jH)) 
should be written as 

^2 (x, y) = Si U.e^'^-^+2^'^' + 5.e-ifc..^-2i0. _ c-Zi-e^"^ - ^e"^"^ J e'^y^ (7) 

where 

Si = sign (V^ — E) ; hki^ = y 2m \E — Vi\ coscpi; hkiy = y 2m \E — Vi\ sincpi 



K.ix = \ Ikf^ + 2kfy] hi= \ \ll + sin^ 0i - sm ( 



To find the transmission coefficient through barrier (j2)), one should set rfi = for x < 0, 
63 = C3 = for X > D and satisfy the continuity conditions for both components of 
the wavefunction and their derivatives. For the case of an electron beam that is incident 
normally {(f) = 0) and low barriers Vq < E (over-barrier transmission), we obtain ipi = —ip2 
both outside and inside the barrier, and the chirality of fermions in bilayer graphene does 
not manifest itself. In this case, scattering at the barrier (jSJ is the same as for electrons 
described by the Schrodinger equation. However, for any finite (even in the case Vq < -E), 
waves localized at the barrier interfaces are essential to satisfy the boundary conditions. 
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FIG. 3: Transmission probability T for normally incident electrons in single- and bi- layer graphene 
(red and blue curves, respectively) and in a non-chiral zero-gap semiconductor (green curve) as 
a function of width D of the tunnel barrier. Concentrations of charge carriers are chosen as 
n = 0.5 X 10^^ cm^^ and p = 1 x 10^^ cm~^ outside and inside the barrier, respectively, for all 
three cases. This yields barrier's height of ~ 450 meV for graphene and ~ 240 meV for the other 
two materials. Note that the transmission probability for bilayer graphene decays exponentially 
with the barrier width, even though there are plenty of electronic states inside the barrier. 



The most intriguing behavior is found for Vq > E, where electrons outside the barrier 
transform into holes inside it, or vice versa. Examples of the angular dependence of T in 
bilayer graphene are plotted in Fig. |2b- They show a dramatic difference as compared with 
the case of massless Dirac fermions. There are again pronounced transmission resonances at 
some incident angles, where T approaches unity. However, instead of the perfect transmission 
found for normally-incident Dirac fermions (see Fig. |2t), our numerical analysis has yielded 
the opposite effect: Massive chiral fermions are always perfectly reflected for angles close to 

= 0. 

Accordingly, we have analyzed this case in more detail and found the following analytical 
solution for the transmission coefficient t: 

4ikik2 



t 



(^2 + ikif e-^^D - {k2 - ikif e^^^ '' 



(8) 



where subscripts 1,2 label the regions outside and inside the barrier, respectively. Par- 
ticularly interesting is the case of a potential step, which corresponds to a single p — n 



junction. Eq (jS)) shows that such a junction should completely reflect a normally-incident 
beam (T = 0). This is highly unusual because the continuum of electronic states at the 
other side of the step is normally expected to allow some tunneling. Furthermore, for a 
single p — n junction with Vq ^ E, the following analytical solution for any (f) has been 
found: 

T=^stn\2<j)) (9) 

which again yields T = for = 0. This behavior is in obvious contrast to single- layer 
graphene, where normally-incident electrons are always perfectly transmitted. 

The perfect reflection (instead of the perfect transmission) can be viewed as another 
incarnation of the Klein paradox, because the effect is again due to the charge-conjugation 
symmetry (fermions in single- and bi-layer graphene exhibit chiralities that resemble those 
associated with spin 1/2 and 1, respectively)^'^. For single- layer graphene, an electron 
wavefunction at the barrier interface matches perfectly the corresponding wavefunction for 
a hole with the same direction of pseudospin (see Fig. la), yielding T = 1. In contrast, for 
bilayer graphene, the charge conjugation requires a propagating electron with wavevector k 
to transform into a hole with wavevector ik (rather than —k), which is an evanescent wave 
inside a barrier. 

COMPARISON WITH TUNNELING OF NONCHIRAL PARTICLES 
For completeness, we compare the obtained results with the case of normal electrons. If 
a tunnel barrier contains no electronic states, the difference is obvious: the transmission 
probability in this case is well known to decay exponentially with increasing barrier's width 
and height^ so that the tunnel barriers discussed above would reflect electrons completely. 
However, both graphene systems are gapless, and it is more appropriate to compare them 
with gapless semiconductors having nonchiral charge carriers (such a situation can be realized 
in certain heterostructure o^?'^^ ) . In this case, one finds 

t = , ^^ 2 , (10) 

where k^ and q^ are x-components of the wave vector outside and inside the barrier, re- 
spectively. Again, similarly to the case of single- and bi-layer graphene, there are resonance 
conditions qxD = nN, N = 0, ±1,... at which the barrier is transparent. For the case of 
normal incidence (0 = 0) the tunneling coefficient is then an oscillating function of tun- 
neling parameters and can exhibit any value from to 1 (see Fig. 3). This is in contrast 
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FIG. 4: The chiral nature of quasiparticles in graphene strongly affects its transport properties, (a) 
A diffusive conductor of a size smaller than the phase-coherence length is connected to two parallel 
one-dimensional leads. For normal electrons, transmission probability T through such a system 
depends strongly on the distribution of scatterers. In contrast, for massless Dirac fermions, T is 
always equal to unity due to the additional memory about the initial direction of pseudospin (see 
text), (b) Schematic diagram of one of the possible tunneling experiments in graphene. Graphene 
(light-blue) has two local gates (dark-blue) that create potential barriers of a variable height. The 
voltage drop across the barriers is measured by using potential contacts shown in orange. 

to graphene, where T is always 1, and bilayer graphene, where T = for sufficiently wide 
barriers D > X. This makes it clear that the drastic difference between the three cases is 
essentially due to different chiralities or pseudospins of the quasiparticles involved rather 
than any other feature of their energy spectra. 

IMPLICATIONS FOR EXPERIMENT 

The found tunneling anomalies in the two graphene systems are expected to play an 
important role in their transport properties, especially in the regime of low carrier con- 
centrations where disorder induces significant potential barriers and the systems are likely 
to split into a random distribution of p-n junctions. In conventional 2D systems, strong 
enough disorder results in electronic states that are separated by barriers with exponentially 
small transparenc y^^ 1^^ . This is known to lead to the Anderson localization. In contrast, 
in both graphene materials all potential barriers are relatively transparent {T ^ 1 at least 
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for some angles) which does not allow charge carriers to be confined by potential barriers 
that are smooth on atomic scale. Therefore, different electron and hole "puddles" induced 
by disorder are not isolated but effectively percolate, thereby suppressing localization. This 
consideration can be important for the understanding of the minimal conductivity ~ e^ /h 
observed experimentally in both single-layer^^ and bilayer— graphene. 

To elucidate further the dramatic difference between quantum transport of Dirac fermions 
in graphene and normal 2D electrons. Fig. E^ suggests a gedanken experiment where a diffu- 
sive conductor is attached to ballistic one- dimensional leads, as in the Landauer formalism. 
For conventional 2D systems, transmission and reflection coefficients through such a con- 
ductor are sensitive to detailed distribution of impurities and a shift of a single impurity 
by a distance of the order of A can completely change the coefficients^^. In contrast, the 
conservation of pseudospin in graphene strictly forbids backscattering and makes the disor- 
dered region in Fig. |3^ always completely transparent, independent of disorder (as long as 
it is smooth on the scale of the lattice constanlji^). This extension of the Klein problem to 
the case of a random scalar potential has been proven by using the Lippmann-Schwinger 
equation (see the Supplementary Information). Unfortunately, this particular experiment 
is probably impossible to realize in practice because scattering at graphene's edges does 
not conserve the pseudospi n^^i^^ . Nevertheless, the above consideration shows that impu- 
rity scattering in the bulk of graphene should be suppressed as compared to the normal 
conductors. 

The above analysis shows that the Klein paradox and associated relativistic-like phenom- 
ena can be tested experimentally using graphene devices. The basic principle behind such 
experiments would be to employ local gates and collimators similar to those used in electron 
optics in 2D gases^^*^. One possible experimental setup is shown schematically in Fig. EJd. 
Here, local gates simply cross the whole graphene sample at different angles (for example, 
90° and 45°). Intrinsic concentrations of charge carriers are usually low (~ 10^^ cm~^), 
whereas concentrations up to 1 x 10^^ cm~^ can be induced under the gated regions by the 
bipolar electric field effect^. This allows potential barriers with heights up to Vq ~ 0.4 eV 
and ~ 0.23 eV for single- and double-layer samples, respectively. By measuring the volt- 
age drop across the barriers as a function of applied gate voltage, one can analyze their 
transparency for different Vq. Our results in Fig. |21 show that for graphene the 90° barrier 
should exhibit low resistance and no significant changes in it with changing gate voltage. In 
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comparison, the 45° barrier is expected to have much higher resistance and show a number 
of tunnehng resonances as a function of gate voltage. The situation should be qualitatively 
different for bilayer graphene, where local barriers should result in a high resistance for the 
perpendicular barrier and pronounced resonances for the 45° barrier. 

Furthermore, the fact that a barrier (or even a single p — n junction) incorporated in a 
bilayer graphene device should lead to exponentially small tunneling current can be exploited 
in developing graphene-based field effect transistors (FET). Such transistors are particularly 
tempting because of their high mobility and ballistic transport at submicron distancea^ i^^i^^ . 
However, the fundamental problem along this route is that the conducting channel in single- 
layer graphene cannot be pinched off (because of the minimal conductivity), which severely 
limits achievable on-ofF ratios for such FETsS. and, therefore, the scope for their applications. 
A bilayer FET with a local gate inverting the sign of charge carriers should yield much higher 
on-off ratios. 

OUTLOOK 

We have shown that the recently found two carbon allotropes provide an effective medium 
for mimicking relativistic quantum effects. On the one hand, this allows one to set up 
such exotic experiments as the one described by the Klein paradox and could be useful for 
analysis of other relevant QED problems. On the other hand, our work also shows that 
the known QED problems and their solutions can be applied to graphene to achieve better 
understanding of transport properties of this unique material that is interesting from the 
view point of both fundamental physics and applications. 
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